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Two-pion correlation functions are analyzed at mid-rapidity for three systems 
(14.6 A-GeV/c Si-I-Al, Si-I-Au, and 11.6 A-GeV/c Au-|-Au), seven distinct centrality condi- 
tions, and different kr bins in the range 0.1-0.5 GeV/c. Source reference frames are determined 
from fits to the Yano-Koonin source parameterization. Bertsch-Pratt radius parameters are shown 
to scale linearly with both number of projectile and total participants as obtained from a Glauber 
model calculation. A finite emission duration that increases linearly with system/centrality is also 
reported. The my dependence of the Bertsch-Pratt radii is measured for the central Si-fAu and 
Au-|-Au systems. The system/centrality dependence is investigated separately for both high and 
low TTiT regions. 
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Bose-Einstein correlations of identical charged pions 
were first observed and used to extract a source size for 
pion emission in pp annihilations j^]. This technique 
has since been applied to heavy ion collisions ranging 
m energy from 1.8 A-GeV/c at the Bevalac [|[ || to 
y^snn — 130 GeV at RHIC |^ (see also recent re- 
views 1^, 1^ and references therein). Symmetrization of 
the two-pion wave function results in an enhancement of 
the two-particle correlation in a region of low relative mo- 
mentum, the extent of which is inversely proportional to 
the size or radius of the emitting source. The technique 
is most commonly referred to as HBT, after the simi- 
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lar technique pioneered by Hanbury-Brown and Twiss to 
measure stellar radii from intensity interferometry with 
radio waves Q . The form of the two-pion correlation de- 
pends upon the emission function assumed for the pion 
source. For static source distributions with no final state 
interactions, such as the Fermi Statistical Model first 
used by Goldhaber et al. Q|, the correlation function is 
related to the square of the 3-dimensional Fourier trans- 
form of the source distribution, /o(r), with respect to the 
pion relative momentum, q = Pi — P21 



C(Pl,P2) 



^'(Pl,P2) 



l + |p(q)|^ 



(1.1) 



^(Pl)^(P2) 

For a multi-dimensional, Gaussian parameterization of 
the source, the enhancement in the relative momentum 
correlation is given by a multi-dimensional Gaussian with 
each Gaussian width inversely proportional to the canon- 
ically conjugate dimension, or radius. For dynamic, 
rapidly expanding sources the radii decrease with increas- 
ing transverse mass, tot, of the pion pair j|, ^ [ill, |l^, 

[Tst , where tot — fcr^ + "t-^, and kx is the mean trans- 
verse momentum of the pion pair. These TOr-dependent 
radii correspond to the relative separations of the pious 
with low relative momentum and are commonly referred 
to as "lengths of homogeneity" . 

Interest in measuring source sizes in heavy ion colli- 
sions comes partly from the expectation that the presence 
of a QCD phase transition will lead to large source sizes 
and long lifetimes for particle emission ^ . The 

absence of such a signal, coupled with the complex and 
subtle nature of the mx dependence, has left many puz- 
zled over how best to interpret the available HBT results. 
However, in heavy ion collisions, there exists additional 
geometric information from the size of the system, the 
centrality of the collision, and also the reaction plane. 
Results on the pion source shape relative to the reac- 
tion plane for non-central collisions |l^ provide one 
important confirmation that HBT radii arc strongly cor- 
related with physical dimensions of the nuclear overlap 
region. A previous comparison of HBT sources to nu- 
clear geometry Jl8[ also confirmed this correlation, but 
the compilation included different source parameteriza- 
tions, online triggers, and experimental acceptances. 

We present a study of the dependence of the radii 
on system size and centrality for the two-pion cor- 
relation function from 14.6 A-GeV/c Si-|-Al, Si+Au, 
and 11.6 A-GeV/c Au-hAu collisions at the BNL AGS 
measured by the E802 rotating magnetic spectrometer 
(Henry Higgins). The radii are compared to three ge- 
ometric quantities: the effective nuclear radius for the 
number of projectile participants, the radius for the to- 
tal participants, and the transverse distribution of binary 
collisions calculated with a Glauber model. The data sets 
of the three systems span a range of seven distinct nu- 
clear geometries. The tot dependence is examined for 
the high statistics central Si-|-Au and Au-|-Au systems 
and compared to measurements of Au-|-Au central colli- 
sions by the E866 Forward Spectrometer. Furthermore, 



the centrality dependence is investigated separately for 
high and low values of mx for all systems. 

Sec. H describes the experimental apparatus relevant 
to this analysis. Sec. Ill covers the data reduction. The 
correlation analysis procedure is described in Sec. ^ and 
results are presented in Sec. ^ The conclusions are given 
in Sec. Vl. Appendix ^ contains the complete set of ta- 
bles for the correlation radii. A Monte Carlo study to 
check systematic effects in the Yano-Koonin fits is pre- 
sented in Appendix Details and systematic studies of 
the Coulomb correction are provided in Appendix M. 



II. EXPERIMENT 

The Si+Al/Au data were collected in 1991 and 1992 
by BNL E859, and the Au-|-Au data were collected in 
1992 and 1994 by BNL E866. Both experiments were ex- 
tensions of BNL E802. This series of experiments ran in 
the Bl line of the Brookhaven AGS from 1986 through 
1995. The complete experimental setups are described 
elsewhere [|l9|, The components essential for this 

analysis are a Target Multiplicity Array (TMA), a Zero- 
degree Calorimeter (ZCAL) a series of four drift 
chambers (T1-T4), two trigger chambers (TR1,TR2), 
two multi-wire proportional chambers (TRF1,TRF2) a 
time of flight wall (TOF), beam counters (BC), and a 
level-2 trigger (LVL2). 

The TMA measured total charged particle multiplicity 
from resistive tube pads configured in a wall and barrel 
array surrounding the target. A sum of the discriminated 
pad signals was used in the level- 1 trigger to enhance both 
central and peripheral data sets. The TMA was present 
only for the Si beam and was removed prior to the Au 
beam commissioning in 1992. The ZCAL consisted of 
alternating layers of Pb and scintillator and was placed 
11 m downstream from the target. It was completely 
rebuilt and incorporated into the level- 1 trigger prior to 
the start of Au beam in 1992. For the data sets presented 
here the rms energy resolution ranged from 2.0y/E for the 
Si beam to for the Au beam. 

Beam definitions and to for TOF were provided by up- 
stream beam counters (BC) in anti-coincidence with a 
series of veto paddles. The interaction trigger condition 
was established by a buUseye (BE) scintillator detector, 
placed just before the ZCAL. For the Au beam, this was 
replaced by a Cercnkov radiator. The interaction trigger 
was defined by Z < 26.5 for the Si Beam, and Z < 73 for 
the Au beam. 

The four drift chambers were placed two before and 
two after the magnet. Tl, T2, T3, and T4, were posi- 
tioned approximately 1, 1.5, 4, and 4.5 m downstream 
from the target, respectively. The drift chambers pro- 
vided a position resolution of 150 /zm, and a minimum 
two-track separation of 2 mm. For this analysis, the mag- 
net was run at 0.4 Tesla, for an integrated field of 0.585 T- 
m. The trigger chambers, TRl and TR2, were located 
behind T3 and T4, respectively. For the 1994 running. 
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two highly segmented muhi-wire chambers, TRFl and 
TRF2 were added behind Tl and T2, respectively, to 
improve track reconstruction in the higher multiplicity 
environment. Resolutions for the TRF chambers were 
approximately 200 fim. The TOF was placed about 6 m 
downstream and particle flight times were measured with 
an rms resolution of 120 ps. 

The level-2 trigger (LVL2) was an integral part of the 
data taking and was especially important for collecting 
two-pion events in peripheral collisions. It consisted of 
TRl, TR2, TOF, and a series of LeCroy CAMAC mod- 
ules used for fast readout, hit storage, memory lookup, 
and logic |Q. The trigger algorithm looped over all 
TRl-TOF hit combinations and matched hits on TR2 
for tracks projected from the target. Momentum and 1/(3 
for valid combinations were determined through lookup 
tables based on hit positions and time-of- flight. One or 
two independent regions in mass vs. momentum could be 
specified in any logical combination to form the trigger 
condition. All events for which the LVL2 did not produce 
a veto within 40 /isec were accepted. A set of events in 
which the LVL2 decision was recorded but not used to 
veto were collected for each LVL2 condition used in the 
experiment. All data sets presented here made use of the 
level-2 trigger, sometimes in combination with a level- 1 
centrality trigger in the TMA or ZCAL. LVL2 rejection 
factors ranged from 3-10 for Au-|-Au central to Si-|-Al 
peripheral events. Trigger inefficiencies were determined 
to be of order 1% for these data, and were dominated 
by the wire chamber inefficiencies. No detectable trigger 
bias was found for the two-pion data sets. 

Nine data sets were taken covering seven different sys- 
tem and centrality conditions. These data sets are sum- 
marized in order of increasing system size and centrality 
in Table ^ The online cut column refers to the multiplic- 
ity hardware centrality definition for the Si beam data 
and the forward energy definition for the Au beam data. 
Sets 1 and 2 are for 7r+ pairs from a 14.6 A-GeV/c Si 
beam incident on a 1.63 g/cm^ Al target (6% beam in- 
teraction length). Set 1 combines both peripheral TMA 
triggered data and minimum bias data. Set 2 includes 
only the minimum bias data. Sets 3-6 are for pion pairs 
from a Si beam incident on a 0.944 g/cm^ Au target 
(1%). Two different TMA level- 1 thresholds were used 
to collect negative pion pairs for peripheral and semi- 
peripheral events (sets 3 and 4). The central TMA trig- 
ger was used to collect pion pairs of both signs (sets 5 
and 6). Negative pion pairs from 11.6 A-GeV/c Au-|-Au 
collisions were taken using the same Au target (1.5% in- 
teraction length for Au beam) for central triggers in 1992 
(set 9), and for minimum bias events in 1994 (sets 7 and 
8). The Au-I-Au minimum bias data set was divided at 
EzcAL= 550 in order to match the online centrality con- 
dition for data set 9. All data sets were subject to offline 
cuts to reject beam pile-up and to improve centrality def- 
inition. Table || gives the final pion pair statistics after 
all offline cuts. 

The final TMA multiplicity and ZCAL forward energy 



TABLE I: Summary of correlation data sets. 









Online 


Offline 




Set 


System 




Centrality 


Centrality 


Pairs 


1. 


Si-I-Al^ 




periph./m.b. 


Ntma<45 


83602 


2. 


Si+Al^ 


27r+-|-X 


min. bias 


Ntma>30 


78713 


3. 


Si-fAm 


27r--fX 


periph. 


Ntma<75 


50713 


4. 


Si-hAm 


27r-+X 


semi-periph. 


Ntma<115 


98468 


5. 


Si-HAu^ 


27r-+X 


central 


Ntma>75 


232296 


6. 


Si-I-Au^ 


27r++X 


central 


Ntma>75 


76661 


7. 


Au+Au- 


27r-+X 


min. bias 


EzcAL>550 


77837 


8. 


Au+Au- 


27r-+X 


min. bias 


EzcAL<550 


88198 


9. 


Au-f Au- 


27r"+X 


central 


none 


85573 



distributions for all nine data sets are shown in Fig. 
Only the ZCAL distributions are shown for the Au-|-Au 
system (sets 7-9) — a new multiplicity array (NMA) 
was in place for the collection of data in sets 7 and 8 
but was not used in this analysis. The negative values 
in some of the Ezcal distributions were the results of 
excessive noise in the ZCAL electronics readout. This ef- 
fect is taken into account in the centrality determination 
described in the next section. 
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FIG. 1: TMA multiplicity and ZCAL forward energy dis- 
tributions for all pion pair events. Central distributions are 
densely hatched, peripheral distributions are unfilled, and the 
Si+Au semi-peripheral distribution is sparsely hatched. Data 
set numbers from Table | are printed alongside the distribu- 
tions. All distributions are normalized to unit area. 
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III. DATA REDUCTION 

A. Centrality Determination 

The forward energy distributions from ZCAL were 
used to determine the collision geometry for each system. 
The number of participant nucleons from the projectile 
is nominally calculated from the ZCAL forward energy 
(Ezcal) and the kinetic energy per nucleon (Enucieon) 
according to, 



Impact Parameter 



Total Participants 



zc 



-^proj — A — 



ZCAL 
lucleoii 



(3.1) 



negative Ezcal val- 
The mean values for 



In applying Eq. 3.1 to the data, 
ues were included in the average. 

^"^^proj^^^ are given in Table ||. 

To account for the effect of the excess noise in the 
ZCAL, and to obtain a more accurate determination of 
the collision geometry, a Monte Carlo Glauber model was 
used to simulate the forward energy distributions. We 
use a value of 30 mb for the N-N cross-section and the 
nuclear distributions are given by a Wood-Saxon distri- 
bution with nuclear radius, r ~ l.OTA^/"^ fm, and surface 
thickness 0.55 fm p^ . 

The ZCAL value for each Monte Carlo event in a given 
system was calculated by sampling from a Poisson distri- 
bution set by the beam distribution for the data. This 
distribution was then smeared with a Gaussian to sim- 
ulate the additional noise observed in the distributions. 
The width of the Gaussian term (10-20 GeV) was chosen 
to match the noise in the target-out corrected minimum 
bias distribution while maintaining agreement for beam 
events. 

The model was used to calculate distributions for the 
number of projectile, target, and total participants and 
the transverse rms distribution of binary collisions (prms) ■ 
Each Monte Carlo event was assigned a weight accord- 
ing to it's calculated ZCAL value, defined by the ZCAL 
distribution for 27r events divided by the minimum bias 
distribution. This weighting procedure is necessary to 
reproduce the correct ZCAL distribution for events with 
two pions in the spectrometer. The weighted distribu- 
tions for the irnpact parameter and total participants are 
shown in Fig. |. The mean values are listed in Table 0, 



The average value of iVproj calculated with Eq. 3.1 



are similar to the values of A^proj calculated using the 
Glauber model. All geometric quantities are observed to 
increase with increasing centrality within each system. 
The geometric values are identical for the two central 
Si-|-Au data sets and nearly equivalent for the two cen- 
tral Au-I-Au data sets. 



B. Pion identification and pair selection 

The tracking algorithm for the 1991-1992 data (sets 1- 
6, 9) |2l began with a road-finder to connect TOF, TRl 




5 10 

b(fm) 



total 



FIG. 2: Distributions for impact parameter (left) and to- 
tal participants (right) from Monte Carlo Glauber model for 
all data sets. More central distributions are more densely 
hatched, as in Fig. |l| 



TABLE II: Average geometric quantities calculated from 
ZCAL distributions and a Monte Carlo Glauber model. Par- 
ticipant numbers are raised to the one-third power and then 
averaged. The value for prms is the transverse rms distribution 
of all nucleon-nucleon interactions (in units of fm) averaged 
over all collisions. 



System 


Set 


proj 


^1/3 
pro.i 


^1/3 
' targ 


^1/3 
^ 'total 


prms 


Si+Al 


1 


L94 


1.76 


1.75 


2.20 


1.15 


2 


2.46 


2.23 


2.18 


2.78 


1.57 




3 


2.19 


2.10 


2.43 


2.87 


1.59 


Si+Au 


4 


2.49 


2.43 


2.94 


3.42 


1.85 


5 


2.87 


2.88 


3.79 


4.29 


2.20 




6 


2.87 


2.88 


3.79 


4.29 


2.21 




7 


4.40 


3.95 


3.97 


5.00 


2.52 


Au+Au 


8 


5.52 


5.43 


5.44 


6.84 


3.35 




9 


5.45 


5.48 


5.46 


6.89 


3.36 



and TR2 hits. Additional hits collected on the down- 
stream tracking chambers were used to construct the 
downstream track segment. This was projected through 
the magnet using an effective dipole field approximation 
assuming a target vertex and hist were collected on Tl 
and T2. All hits associated with the track were fit with 
the vertex constraint removed. Tracks were required to 
pass goodness of fit cuts and to project to within 2 cm of 
the nominal target position. The single particle momen- 
tum resolution was determined by GEANT simulation 
and parameterized as Cp = ap/l3 (B bp^ , with a = 1.2% 
and b = 0.6% (GeV/c)"i. For the 1994 data (sets 7- 
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8) upstream track segments were found and fit indepen- 
dently using the additional hits on TRFl and TRF2. 
The two track segments were then fit using the same 
magnetic field approximation, and applying the same 
vertex and goodness of fit cuts. Using additional in- 
formation from TRFl and TRF2 the momentum reso- 
lution for 1994 was improved to give a = 0.85% and 
6 = 0.12% (GeV/c)-i [11 . 

Pions were identified up to a momentum of 1.85 GeV/c 
by requiring a unique association with a TOF hit, and a 
measured 1//3 that lies within 3a of the expected value 
for pions and outside of the 3a band for kaons. The 
dominant background comes from electron contamina- 
tion, determined to be less than 5% |2^ in the region 
0.54 < p < 1.3 GeV/c. Pion pair momentum distribu- 
tions for the central Si-|-Au and Au-|-Au data sets are 
shown in Fig. ^. The Si-beam data were taken with a 
spectrometer setting of 14 degrees, which covers an an- 
gular region from 14 to 28 degrees and leads to a mean 
pair-rapidity for pions approximately equal to Fnn=1-72. 
For the Au-beam the spectrometer was rotated to 21 de- 
grees, primarily to achieve lower track densities in the 
upstream chambers. The mean pion pair-rapidity for the 
Au-beam data is 1.43, slightly backwards of Inn = 1.61. 



parameterized by a two-dimensional Gaussian, 



= 1 -0.86exp 



[2.75)^601 + 561 



2(0.0121)2 



PIxIPIz- P2x/P2z, 
Ply/Plz -P2y/P2z- 



(3.2) 



A detailed GEANT simulation of two-pion events pro- 
duced the same Gaussian widths for the two-track effi- 
ciency, but the coefficient for the exponential was closer 
to unity. One expects this coefficient to be less than 
unity in the data due to the presence of mi s- as sociated 
hits from other tracks. The inverse of Eq. 3.2 was ap- 
plied as a correction in the region where the efficiency 
was greater than ~50%, defined as. 



^{2.75)^601 + 60. 



I > 0.011. 



(3.3) 



Pairs not in this region were removed from the analysis. 

For the 1994 data (sets 7 and 8), the two-track effi- 
ciency is flat but then drops sharply for close tracks. For 
these data, no correction was applied and the two-track 
cut is given by, 



691 



91 > 0.008. 



(3.4) 
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FIG. 3: Distributions in kr and rapidity for Si+X (14° spec- 
trometer setting) and Au-|-Au (21° setting). Arrows in- 
dicate rapidity for nucleon-nucleon center of mass for 
14.6 A-GeV/c Si and 11.6 A-GeV/c An beams. 



The two-track resolution was determined separately for 
the 1991-1992 and 1994 tracking algorithms, by examin- 
ing the ratio of signal to event-mixed pairs versus the 
relative slopes 69 x and 69y upstream of the magnet. For 
the 1991-1992 data set, the two-track efficiency is well 



IV. CORRELATION ANALYSIS 

A. Parameterizations 

Ideally the Bose-Einstein correlation function is de- 
fined to be the ratio of the probability for emitting two 
pions in the presence of Bose-Einstein symmetrization 
to the probability for the case with no symmetrization. 
The latter can be approximated using Monte Garlo or 
event mixing ||2^ techniques. For each data set, an event- 
mixed background was formed from pion pairs from dif- 
ferent events chosen at random from the entire data set. 
The number of background pairs was chosen to be 5-10 
times the signal to achieve statistical errors smaller than 
those for the signal, but still Poisson distributed [ p6| . 
This background differs from the ideal in two respects, 
it contains a residual correlation distortion first noted by 
Zajc and it lacks two-body final state interactions 
such as the Coulomb interaction. The effect of the resid- 
ual correlation is included in the systematic error esti- 
mates. We correct for the Coulomb interaction using an 
analytic approximation for the Coulomb wave function 
of two-pions emanating from a finite static source. This 
source is assumed to be a spherically symmetric Gaussian 
in the pair center-of-mass frame. This approximation has 
been described before |2j| and is reviewed in Appendix H. 



The two-track cuts (Eq. 3.3 and 3.4) are applied to both 
the signal and background. Other final state interactions 
such as the di-pion strong interaction p9|, 30 and the 
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Coulomb interaction with the nuclear remnant are ex- 
pected to be negligible |^ and are not considered. The 
correlation function is then given by, 



1 



1 



(4.1) 

Here, A and B are the signal and event-mixed back- 



ground, e is the two-track efficiency defined by Eq. 3.2, 
and F is the Coulomb correction defined by Eq. Bl. 
Qinv is the magnitude of two-pion relative 4-momentum 
in the pair rest frame, and Rinv is the parameter returned 
by the 1-D invariant Gaussian parameterization. 



(4.2) 



Fits to this parameterization are used to perform an it- 
erative determination of the Coulomb correction. Begin- 
ning with the assumption of a point source (i.e. Gamow 
correction), the i?inv parameter converges to within 1% 
after two iterations for all data sets. 

The Yano-Koonin parameterization |^ is used to ex- 
tract the longitudinal velocity of the source, 

C2(Qi„v,gi) = l + Ae-(«^+-^)^=^(^«-'^='')^+('3"^vH^), 

(4.3) 

where qi = Pzi — Pz2 and (3s is the longitudinal compo- 
nent of the source velocity. R and r are the spherically 
symmetric Gaussian radius and emission duration, re- 
spectively. 

Using the source frames determined from the Yano- 
Koonin fits, all data sets were fit to the 3-D Bertsch-Pratt 
parameterization |l^ , 

( 2 r>2 2 n2 2 n2\ 

C2{quqs,qo) = 1 + Ael-''' ^' -«ofl„) , (4.4) 

where q\ is the long itud inal component of the relative mo- 
mentum from Eq. 4.3. The transverse relative momen- 



tum is separated into components perpendicular (qs) and 
parallel (qo) to the pion pair velocity. This parameteriza- 
tion measures the longitudinal and transverse correlation 
lengths, R\ and i?s, respectively. For an azimuthally sym- 
metric source without dynamical correlations, the pion 
emission duration is calculated from i?o, ^s, and the 
transverse velocity of the pion pair, (3t, 



T = ^jRo'-Rsy <Pt > 



(4.5) 



A variation of this parameterization suggested by Chap- 
man [Q, includes a cross-term, Rio, which is expected 
to be non-vanishing for sources that are asymmetric in 

z, 

C2{quqs,qo) = 1 + Ael-''?^?-^^^^ -'o«o-2,,.ofiy _ (4.6) 

It is instructive to relate the cross-term to the Gener- 
alized Yano-Koonin parameterization (GYK) using 
the relation, qo — q- (3pr. 

C2{quqT,qo) = 1 + Ae(-'?^?-«T«^-9^^) (4.7) 
= 1 + Ae(~'?'(^'+ft'^')-«T(^T+/3r^')-2<2i<2oA/3rr"))_ 



This provides a simple explanation for the cross-term, 
Rio in the absence of dynamical correlations, and illus- 
trates why the cross-term vanishes for sources which are 
symmetric about the z = plane. 



B. Minimization 

A new log-likelihood minimization function is used in 
this analysis. It assumes Poisson probability distribu- 
tions for both signal and background, with means given 
by /i and v respectively. The correlation function is de- 
fined as the ratio of these means by imposing a delta- 
function constraint on the conditional probability for C, 
given A and B, 



P{C\A,B) = / dfidiy 



A\ 

{A + B^ 



B\ 



D! 



A\B\ 



(C+1) 



A+B+2 ■ 



(4.8) 



Taking the negative of twice the log and collecting lead- 
ing order terms in A and B leads to the minimization 
function used in this analysis. 



A\n 



C{A + B) 

A{C + l) 



Bin 



A + B 
B{C+1) 



Note that Eq. 



(4.9) 



approaches a chi-squared distribution, 
■^— — in the limit of large A and B. We compare 
Eq. [4.9| to another log-likelihood minimization function 
commonly used IsJ, 



2 

XpML 



-2{C2B 



'C2B 



)ln 



{C2B -A) + [a 



(4.10) 



'C2B/ 



'C2B) 



This form also assumes a Poisson distributed signal for A, 
but the error terms for the background have been derived 
by working backwards from the chi-squared limit. All 
terms in Eq. |4.9| follow naturally from the assumption of 
Poisson distributions for A and B and the log-likelihood 
method. For this reason, the new minimization function 
may be more accurate for fi ts in volv ing b ins with small 
counts. However, both Eqs. 4.9 and 4.11 were shown to 
give consistent results in fits to the data. The minimiza- 
tion was performed using the MINUIT package |33] . 



V. RESULTS 

For each of the nine data sets in Table | the three pa- 
rameterizations, Qinv, Yano-Koonin, and Bertsch-Pratt 
were performed in order. The Qinv fits were used to de- 
termine the Coulomb correction. The longitudinal veloc- 
ity of the source frames were obtained from the Yano- 
Koonin fits, and the HBT radii from the Bertsch-Pratt 
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fits were used in analysis of centrality, system, and kT de- 
pendence. A complete set of tables for all fitted radii are 
given in Appendix |^. 



fits are summa- 



A. Qinv Fit Results 

Fit results f or the final it eratio n Qi, 
rized in Tabl e |yill| . Tables |VIIl| through ^ are found 
in Appendix |C[. These radii were taken as input to the 
Coulomb correction calculation for all subsequent multi- 
dimensional parameterizations for these data sets. For 
eac h sys tem, i?inv increases with collision centrality. Ta- 
ble VII] also lists the mean transverse mass,< uit >, for 
each of the nine data sets. In calculating this mean, it is 
common to restrict the average to only those pairs with 
low relative momentum, which carry the information on 
the source dimensions. However, for systems that vary in 
size, the region of low relative momentum will also vary 
inversely with the source size. To account for this effect, 
the average restricted to the region Qinv ^ '^l Rim using 
the step function O, 



< rrirp > = 



(5.1) 

The quantity < > should be sufficiently general to 
be used over a wide variation of system sizes and exper- 



imental acceptances. Table VIII shows that < rrh^ > is 



never less than 84% of < ttit > for these data. The fully 
corrected Qinv correlation functions fits are shown in the 
left set of panels of Fig. |[ 



Yano-Koonin Fit Results 



By construction, the R and r parameters of Eq. 4.3 
are in the longitudinal rest frame of the source, and are 
invariant under longitudinal boosts. However, large lon- 
gitudinal boosts yield q-distributions that are extended 
along a narrow ridge of the q\ — qq axis, and this presents 
a challenge for binning and minimization. This effect was 
studied with a Monte Carlo simulation and it was deter- 
mined that there is a slight systematic bias in the fitted 
T and Ps when the reference frame differs from the true 
source frame by more than a half unit of rapidity. The 
results of this study are given in Appendix ^. 

To minimize this bias, all data sets were initially fit 
in the Fnn frame, 1.72 for the Si-beam and 1.60 for Au. 
Data sets with values of (3s which differed significantly 
from zero were then refit in the frames expected to yield 
/3s = 0. Only the Si-|-Au central data (sets 5 and 6) 
and Au+Au data (sets 7-9) required a second fit. The 
weighted means for /9s for the Si-|-Au central and Au-fAu 
systems were — 0.24 ± 0.03 and —0.15 ± 0.05 respectively. 
This led to the choice of the Y = 1.50 frame from for 
central Si-|-Au, and the F = 1.45 frame for Au-|-Au. Fi- 
nal fit frames and all parameters for the Yano-Koonin 
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FIG. 4: Correlation functions for data sets 1-9. Data and 
parameterized fits for Qinv (*) are shown in left panels. Right 
panels show slices of the correlation functions and fits for 
Bertsch-Pratt parameterization (t) for g < 40 MeV/c in the 
orthogonal variables. 



These frames are used for all 



fits are listed in Table [X. 
remaining fits. 

All final values for Ps are within two standard devia- 
tions of zero. As with the Qinv results, the Yano-Koonin 
radius increases monotonically with centrality for all sys- 
tems and from smaller systems to larger. The emission 
duration shows a similar but weaker dependence on sys- 
tem size and centrality. We note that because the accep- 
tance for the Si-beam pairs centered at Ynn (see Fig. 13) , 
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the final values for /3s for the Si+Al and non-central 
Si+Au systems are consistent with both the assumptions 
of a fixed source at Inn and longitudinal boost invari- 
ance. The final values of /3s for the symmetric Au+Au 
system are closer to the mean pair rapidity of y = 1.43 
than to Ynn and therefore imply a source that is longi- 
tudinally boost invariant. However, the values of /3s for 
the asymmetric Si-I-Au central system are indicative of 
a pion source shifted backwards to Y = 1.50 and are 
inconsistent with the assumption of longitudinal boost- 
invariance. 



C. Bertsch-Pratt Fit Results 

The results of the Bertsch-Pratt fits with and with- 
out the Chapman R\o cross-term are listed in Ta- 
bles ^ and X], respectively. The system and central- 
ity dependence of the radii will be presented graphi- 
cally in Sec. VD. The R\o cross-terms are consistent 



with zero for the symmetric systems measured at mid- 
rapidity, as expected. The largest values are for the 
Si-|-Au and Au-|-Au central systems, and only for the 
central Si-|-Au— > 27r^ does the R\o term differ from zero 
by more than 2a. 

Slices of the correlation functions for the Bertsch-Pratt 
fits (no cross-term) for all systems in Table ^ are shown 
in the right three panels of Fig. ^. The slices in each 
q-component are averaged over qi < 40 MeV/c in the 
orthogonal q-variables. 

From the values of Ro and Rs in the Bertsch-Pratt fits 
we have e xtra cted values for the emission duration, r, 
using Eq. 4.5. These are compared to the correspond- 



ing emission duration parameters measured in the Yano- 
Koonin fits. Fig. |^ shows the two parameterizations give 
similar results, but with Yano-Koonin r is smaller for 
the lighter systems. This trend is consistent with the 
assumption of spherical symmetry in Yano-Koonin, and 
the trend towards greater asymmetry {R\>Rs) in the 
lighter systems. The Generalized Yano-Koonin (GYK) 
suggested in [|l^ , which separates the transverse and lon- 
gitudin al c omponents of the source, is a logical extension 
to Eq. Use of this parameterization in the future 

may lead to an improved correspondence between emis- 
sion durations determined by the two methods. For the 
remaining sections r is extracted from Eq. 4.5 and its er- 



ror incorporates the Ro-Rs covariant term from the error 
matrix. 

As a systematic check on the reference frames, Bertsch- 
Pratt fits were also performed in the Longitudinal Co- 
Moving System (LCMS) frame, defined as the local frame 
in which the longitudinal pair momentum is zero. These 
results are given in Table XII. All radii are consistent 



with the fixed frame fits to within two standard devia- 
tions. Large differences are not expected given that the 
fixed frames used in the previous fits were consistent with 
the pair LCMS for all but the Si-|-Au central system, 
and the rms width of pair rapidity distribution is small 
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FIG. 5: Comparison between r determined from Yano-Koonin 
and Bertsch-Pratt parameterizations. The dashed line indi- 
cates Tbp — Tyk- 



TABLE III: 

■^totai determined from Glauber model calculations. 



Fit 


Ri 


Rs 


Ro 


r 


a 

A''proj b 


1.44±0.21 
0.40±0.07 
18.2/7 


0.63±0.31 
0.63±0.10 
8.4/7 


1.48±0.19 
0.79±0.07 
10.5/7 


2.01±0.23 
0.45±0.10 
5.6/7 


a 

Motal b 

xVdof 


1.27±0.22 
0.34±0.05 
12.6/7 


0.31±0.34 
0.54±0.08 
2.8/7 


1.43±0.19 
0.59±0.05 
6.3/7 


2.05±0.22 
0.32±0.07 
6.3/7 



M.i). 



D. Centrality and System Dependence 



We examine the dependence of the HBT radii on sys- 
tem size and centrality by comparing the Bertsch-Pratt 
radii given in Table ^ to the geometric quantities of Ta- 
ble ||. Fig. H shows each of the radius parameters and 
the extracted emission duration plotted vs. mean values 
The geometric dependences were 
a + bx. The extracted slopes and 
intercepts are given in Table III. Rg is also fit to to a 
linear dependence on prms, yielding a = —0.01 ± 0.38, 
b = 1.19 ± 0.18, with a xVdof of 3.5/7. All fitted pa- 

1/3 

rameters exhibit a linear dependence on both -/Vp^-oj and 



of iVpYo'j and N^il, 
fit to a linear form 



1/3 



total ■ 



Converting Rs vahies to 3D rms radu (VS X i?s) 



yields slopes of 0.99 ±0.14. and 1.20 ±0.17 for N^^^l^ and 



1/3 

^proj ' respectively. The latter slope is consistent with the 
Ap^'^ scaling noted by J. Bartke [I 



where Ap denotes 



9 



4 

2 

\ 

4 
2 






R / - 

^ — \ \ \ \ \ V- 


1 1 1 1 1 h-. 

' R y ' 

ym^A ^ rms 

-J^ • N ■ - 

^ pro] _ 
"total _ 


X i 



(N X or (N X or p 

^ pro/ ^ tot' "r 



Bertsch-Pratt radii and emission duration vs. projec- 
proj)' total participants (N^^^^^), and the 



FIG. 6 

tile participants (N^^^^), total participants (A'^ 
transverse rms distribution of binary collisions (prms)- AH 
radii are fit to a linear dependence of the form a + bx, with 
fit values given in Table III. 



the atomic number of the projectile. This scahng depen- 
dence was taken from an analysis of mean-free paths of 
nuclear projectiles in emulsion 1 36 , 3^ . It is equivalent to 
the A-dependence of a hard sphere nuclear radius, and is 
consistent with measurements of the nuclear charge dis- 

r- — 1 1/3 

tribution in electron scattering [E3l . The iV^Qtai slopes for 
R\ are 60% smaller than for Rs, whereas for Ro they are 
10% larger. In addition, the t parameter also displays a 
significant slope that is consistent with the dependence 
of Ri. 

We note that the rms values in reference ||l8| were de- 
rived from a variety of trigger conditions, parameteri- 
zations and regions of tot, whereas in this analysis the 
parameterizations are uniform, the trigger conditions are 
accounted for, and the value of < > is in the range 
of 300-325 MeV for Si-hAl/Au and 350-360 MeV for 
Au-|-Au. The implication of the higher < > range 
for Au-I-Au will be discussed in the next section. 



1 /3 

The N^^^^^-^ fits have a slightly lower x'^/dof than the 

^proj fits. The success of the ^totln ^pmj and prms fits 
illustrates that there is more than one way to map the 
system/centrality onto geometry, and that the scaling of 
HBT radii need not be restricted to an effective radius of 
an assembly of wounded nucleons. 

To translate the slopes and intercepts into a true ge- 
ometric dependence of the pion source sizes requires a 
careful study of the mx dependence of each radius. How- 
ever, if the Tot dependence is largely independent of sys- 



tem/centrality, one can extrapolate the entire scaling to 
lower Tot. Then the scaling observed at moderate values 
of Tot will also be valid at lower tot value, where (for 
Rs) the geometric source is better approximated. 

The intercepts for all Rs fits are consistent with zero, 
as one would expect for the dependence of a transverse 
overlap region that is extrapolated to zero participants. 
The intercepts for all other radius parameters differ sig- 
nificantly from zero. This is physically meaningful in 
the context of the standard interpretations given to these 
parameters, where Ri is indicative of the hadronization 
time, and Ro and r are sensitive to the emission duration. 
Thus, even when extrapolating to zero participants, one 
may expect to measure finite lengths and durations for 
the pion hadronization. 



E. Transverse Momentum Dependence 

To study the transverse momentum dependence the 
central Si-|-Au^ 2tt~ data set was subdivided into 
three equal-statistics bins in tot (sets a,c,e), and the 
Si+Au^ 27r+ data set was subdivided into two equal 
bins (sets b,d). The two central Au-|-Au^ 27r~ data sets 
were merged, and then divided into three equal bins (sets 
f,g,h). For this merged data set, the two-track cut of 



Eq. 3.4 was increased from 8 mrad to 15 mrad, and no 
two-track correction was applied. Fits to the full merged 
data set were performed for all parameterizations above 
and were found to be consistent with fits to the sepa- 



rate data sets. The data sets are given in Table XIII , 
along with the fit parameters and values for < to^ >. 
The fully corrected Qinv and Bertsch-Pratt correlation 
functions and fits are shown in Fig. ^. 

The transverse momentum dependence of the two-pion 
correlation function has been discussed by many authors, 
§, |38[ m, 1^, H, and in principle can be used 
to test dynamic models for rapidly expanding sources. 
We have found that the small number of kx bins are 
insufficient to constrain these models. Instead we use 
an exponential form, gio^+P^^T) provide an empirical 
description of the tot dependence of the Si+Au system, 
and use this to check the Au-|-Au system for consistency. 
Values for a and (3 are given in Table IV. Fig. || shows 
the Tot dependence for sets a-h, and for the Si-|-Au radii 
rescaled from N^^^l^ = 2.87 to A^jYtli = 6.89 to match 
Au-|-Au. The solid line shows the tot exponential fit 
also rescaled in this way. The Au+Au data points are 
consistent with the rescaled dependence for each radius. 

Also shown are the E866 forward spectrometer results 
for 7r~ pairs from 1.9 < Y < 2.3 for a centrality range 
that corresponds to the top 24% of the interaction cross- 
section gated on multiplicity Q . These data are less 
central than the Au-|-Au data set at mid-rapidity, which 
correspond to the top ~15% of the interaction cross- 
section gated on forward energy. For comparison, the 
forward spectrometer results are also shown rescaled up 
from iV. 



1/3 
total 



6.3, as determined from the Glauber model 
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FIG. 7: Correlation functions for rriT data sets a-h. Data 
and parameterized fits for Qinv (ir) axe shown in left panels. 
Right panels show slices of the correlation functions and fits 
for Bertsch-Pratt parameterization (t) for q < 40 MeV/c in 
the orthogonal variables. 



TABLE IV: Exponential fit parameters for rriT dependence in 
central Si+Au. 



Variable 


a 




xVdof 


Ri 


1.78±0.21 


-2.24±0.61 


1.7/3 


Rs 


1.53±0.27 


-1.47±0.80 


1.1/3 


Ro 


1.85±0.16 


-1.64±0.49 


2.3/3 


r 


1.59±0.64 


-1.76±2.05 


4.9/3 



ward and mid-rapidity. 
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FIG. 8: rriT dependence of Bertsch-Pratt radii for Si+Au, 
Au+Au, Au+Au at forward rapidity |44[, and for Si+Au 

1/3 — 

rescaled by N^^^^^^ to match Au+Au. The solid line is a fit 
of the form e'""*''""^' to the Si+Au points, rescaled according 

1/3 

the Ai'jjjjjjj for Au+Au. The Au+Au forward rapidity data 
are also shown after being rescaled to match the mid-rapidity 
centrality condition. 



calculations for the 24% most central collisions. Unlike 

1/3 

the determination of N^^^.^^ for the other data sets, this 
was not done using the forward energy distributions for 
events containing 2-pions. This effect is largest for less 
central distributions, and therefore the rescaling of the 
forward rapidity radii in Fig. || should be taken as an 
upper bound. Overall, there is good agreement between 
the the rescaled Si+Au parameterized uit dependence 
and the m-r dependence of the Au+Au radii at both for- 



The parameterized < > dependence of Si+Au radii 
can now be used to check for errors arising from the 20% 
higher < > values of the pion pairs for Au+Au used 
in the centrality analysis. The < > dependence is 
used to interpolate all Bertsch-Pratt radii in Table ^ to 
a common value of < > = 0.325 GeV and the linear 

dependence on iV^^tai and N^ioj is refit. The results are 
given in Table V. The slopes for all parameters except r 
increase slightly. All difference are less than the reported 
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TABLE V: Linear fits of Bertsch-Pratt radii to sys- 
tem/centrality dependence interpolated to < > = 
0.325 GeV. 



Fit 


Ri 


Rs 


Ro 


T 


a 

A^'proj b 

xVdof 


1.23±0.21 
0.47±0.07 
U.9/7 


0.45±0.31 
0.69±0.10 
7.0/7 


1.29±0.20 
0.86±0.07 
8.4/7 


1.91±0.23 
0.49±0.10 
7.7/7 


a 

A^'total b 

xVdof 


1.11±0.22 
0.38±0.06 
9.1/7 


0.16±0.34 
0.57±0.08 
2.1/7 


1.31±0.20 
0.62±0.05 
11.9/7 


2.01±0.23 
0.33±0.07 
10.5/7 



uncertainties, and the observations of the preceding sec- 
tion with respect to slopes and intercepts remain vaHd. 



F. Centrality and Momentum 

To investigate the centrality dependence for low and 
high regions of m-r, each of the data sets 1-7, and 8+9 
combined 'was divided into two mx bins of equal statis- 
tics. The < >, Rinv, and Bertsch-Pratt radii for these 



fits are listed in Tables XIV and KV. The system and 



centrality dependence of the high and low niT regions 
are shown in Fig. ^, plotted only for the A^totai scaling. 
Independent linear fits to the two rriT regions are un- 
able to resolve significant differences in the slopes and 
intercepts for the two regions. The data are compared 

1 /3 

to the A^totai linear fits (solid lines) for the full data sets 
that were shown in Fig. pi The dashed lines are equal to 
the solid line fits interpolated up to 415 MeV and down 
to 275 MeV using the parameterized ttit dependence of 
Si-hAu. 

The high (low) tot values for Ri and Rs fall below 
(above) the linear fits to full mx regions, and are ap- 
proximately consistent with the interpolations to higher 
(lower) Tot. For Ro the high tot points are closer to the 
solid line, and the low tot points fall slightly lower. For 
the T parameter, nearly all points fall below the unbinned 
scaling. One possible explanation is a time-ordered pion 
emission, in which the high tot pions are emitted at a 
time that is displaced relative to the emission of the low 
Tot pions. However, the potentially complex interplay 
between space, time, and momentum must be considered 
fully. At the very least, these results suggest a direc- 
tion for future theoretical work and analysis with higher 
statistics measurements. 

The ratio Ro/Rs has been suggested as way to search 
for changes in the pion emission duration |15| i n heavy 
ion collisions. This ratio is plotted in Fig. |10| for the 
unbinned Si+Au central 27r+ (set 5) and Au-|-Au central 
data (set 8-1-9), and for these data sets subdivided into 
three tot bins. The values of i?o/R-S vary from 1.2 to 
1.5, but do not indicate any clear dependence on the 
Tot bin size. It may be interesting to repeat this study 
in the future with larger data samples and' statistically 




FIG. 9: Bertsch-Pratt radii and emission extracted dura- 
tion vs. total participants {N^^tLi) ^'^^ high "It (A) and 
low rriT (a) bins. The solid lines are linear fits to the full 
rriT range, and the dashed lines are the same linear fits 
rescaled to < > of 275 MeV/c and 415 MeV/c. 



independent bins of different sizes. 



G. Systematic Errors 

The systematic error in determining the centrality is 
expected to be dominated by the ZCAL noise. In the 
course of modifying the noise components in our model, 
we found nominal variations in the derived mean quan- 
tities to be consistent with a systematic error of approx- 
imately 5%. 

The sources of systematic errors for the HBT radii are 
summarized in Tabic VI, and nominal values are given. 
The error due to finite momentum resolution was exam- 
ined for the i?inv fit of the Au-I-Au central data set p5[ , 
and the effect of residual correlations in the background 
is reported elsewhere for the i?inv fit to the central Si-f Au 
data set ^ . The variation for the resolution smear- 
ing applies equally to each of the Bertsch-Pratt radii. All 
other variations are from i?inv parameterizations of the 
central Si+Au— > 27r^ data set. The sign of the variation 
indicates the direction of the effect of incorporating the 
correction/change. Although some of the systematic er- 
rors enter with different signs, the different contributions 
are added in quadrature to yield a total systematic error 
of 3% for the radii and 5% for A. 

The effect of a partial Coulomb correction is considered 
separately in Appendix ^ If a significant fraction of pairs 
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FIG. 10: The ratio of Ro/Rs for the high statistics data sets: 
central Si+Au (set 5, circles) and Au+Au (sets 8+9, trian- 
gles). The open symbols are for the full mr range. The filled 
symbols are for the subdivision into three bins for Si+Au (sets 
a, c, e) and Au+Au (sets g, f, h). 



The other parameters have finite intercepts of 1-2 fm for 

^^proj' ^*total 

Yano-Koonin fits were used to determine the longitudi- 
nal source velocities (/3s) for all systems. For all symmet- 
ric and non-central Si+Au systems measured at Y^n , /Jgis 
consistent with both a fixed source at Fnn and a source 
which exhibits longitudinal boost invariance. The /3s for 
the Au+Au systems are only consistent with boost in- 
variance, and for Si+Au central systems they imply a 
source that is shifted backwards of Inn ■ 

The Tot dependence of the radii has been parameter- 
ized for central Si+Au and shown to be consistent with 
the central Au+Au mx dependence at both mid and for- 
ward rapidity. For the first time, the systcm/centrality 
dependence was investigated independently for both high 
and low tot regions. 

The dependence of R\ and Rs with A'totai at high and 
low Tot is consistent with interpolations based on the pa- 
rameterized Tot dependence of the Si+Au central data. 
However, the i?o and r radii are inconsistent with this 
interpolated tot dependence, implying that these two 
parameters are sensitive to the width of the tot bin. Fi- 
nally, the ratio i?o/ Rs is shown to be above unity for both 
the Si+Au and Au+Au central systems, and binning the 
data in tot does not change this result. 
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TABLE VI: Sources and nominal values of systematic errors 
for the HBT radii and A. 



Source 


R 


A 


momentum resolution 


+2% 


+4% 


5 MeV bins 


±1% 


±2% 


residual correlation 


-2% 


+2% 


two-track correction 


±1% 


±2% 


Coulomb iteration 


+0% 


-1% 


Coulomb resolution 


-1% 


-1% 



(> 50%) are not subject to the Coulomb interaction, then 
this error would dominate for both R and A. 



VI. CONCLUSION 

These data demonstrate unambiguously the strong 
connection that exists between the space-time measure- 
ments from Bose-Einstein correlations and collision ge- 
ometry. All radii and lifetimes were observed to depend 
significantly on several measurements of system size and 
centrality. At < > = 0.325 GeV, Rs was observed 
to scale with the effective rms radius for the total num- 
ber participants, as well as the linear transverse size in a 
Glauber model. For the Ri and r parameters the linear 

1/3 

dependence on N^^^^^^ is 60% weaker, and for R^ it is 10% 
stronger. Only Rs extrapolates to an intercept of zero. 
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APPENDIX A: FRAME DEPENDENCE OF THE 
YANO-KOONIN PARAMETERIZATION 



To study systematic effects of the Yano-Koonin param- 
eterization a Monte Carlo sample of pion events was gen- 
erated for a 3 fm spherically symmetric Gaussian source 
with lifetime 2 fm/c a,t Y — 1.75 i n th e lab. These sim- 
ulated pairs were then fit to Eq. 4.3 while varying the 



fit frame (in which the q-variables were calculated) from 
Y — 0.5 to Y = 3.0 in steps of 0.25 units of rapidity. In all 
fits, the bin sizes were fixed at 20 MeV for q\ and 15 MeV 
for Qinv and go- Fig. ^ shows the results of this study 
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for the fitted radius, lifetime, and source frame. There 



I, 




0.25 



1.75 



3.25 




3.25 



0.25 



1.75 



3.25 

/ 

lab 



FIG. 11: Systematic dependence of Yano-Koonin fit param- 
eters on reference frame for Monte Carlo Gaussian sources, 
with _R = 3 fm, r = 2 fm/c and a source velocity correspond- 
ing to y = 1.75. AF refers to the difference between the fit 
rapidity and the rapidity of the Monte Carlo Gaussian source. 

is a systematic shift in the source frame that is approx- 
imately 10% of the difference between the fit frame and 
the actual source frame. In the simulation this differ- 
ence is noticeable for differences as small as 0.25 units in 
rapidity. The change in source frame is correlated with 
a reduction in the measured lifetime, although R is ob- 
served to be stable for the region tested. 



APPENDIX B: COULOMB CORRECTION 

The Coulomb correction function used in this work is 
an analytic approximation to the full Coulomb wave cal- 
culation. The form of this correction is given by, 



G(77)x 



1 + r]P5{s 



1 + 



(Bl) 



where G is the standard Gamow correction for a point 
source, and P5 is a fifth order polynomial in the dimen- 
sionless quantity s. 



2t:'I] 



P5(s) 



g27rr, _ I 
2fia/ Qinv 
Rinv Q inv / ^ 

min(7r, 8(s - 55^9 + 52sV225)) (B5) 



(B2) 

(B3) 
(B4) 



where, fi is the reduced mass of the pion pair and a is 
the fine structure constant. The full form of the integral 
equation that was solved is described elsewhere [|8| , and 
the FORTRAN code for the tabulated function F{s) is 
available upon request ESI. 



The accuracy of this analytic approximation (Eq. Bl ) 
is demonstrated in Fig. ^ by comparing the correction 
function for a 5 fm source to the Coulomb Wave inte- 
gration calculated by the Correlation After Burner code 
(CRAB) 1 49 . Fig. |lj shows a significant difference be- 
tween the full Coulomb correction and the and Gamow 
correction, but the fast approximation shows good agree- 
ment with the calculation from CRAB. 




- Gamow 

- Eq. Bl (5 fin) 
• CRAB (5 fin) 

- binned correction 
smeared correction 



40 



60 



80 



100 



Qinv (MeV) 

FIG. 12: Coulomb corrections for a point source (Gamow), a 
5 fm spherically symmetric Gaussian calculated with (CRAB) 
and Eq. Bl. The dashed histogram shows the effect of av- 



eraging Eq. Bl over the Si-|-Au data (set 5), and the dot- 



dashed histogram incorporates the effect of momentum reso- 
lution smearing. 



The effect of resolution smearing was estimated by 
evaluating Eq. Bl for the event-mixed background, but 
applying the Coulomb weight at the value of Qinv re- 
turned by a spectral response function. The form of 
the response function was parameterized in momentum 
and polar angle spread according to Monte Carlo single 
particle distributions |Q. Fig. |l^ shows the Coulomb 
correction histogram — weighted background divided by 
background — with and without the response function. 
The effect of the response function is evident only in the 
lowest bin (0-10 MeV/c), and introduced a 1% reduc- 
tion in the 3D HBT radii and A values. Given the small 
size of the effect, it was incorporated into the overall sys- 
tematic error estimate, but the correlation fits used the 
unsmeared correction. 

If the deviation of A from unity is caused by a sub- 
stantial fraction of pions from long lived resonances, it is 
possible that the effect of the Coulomb interaction will 
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TABLE VII: Bertsch-Pratt fit parameters for full, half, and 
no Coulomb correction. 



Parameter 


full 


half 


none 


A 


0.645±0.045 


0.397±0.019 


0.354±0.020 


Ri 


2.96±0.23 


2.53±0.15 


2.32±0.17 


Rs 


2.30±0.24 


2.21±0.21 


1.93±0.23 


Ro 


3.34±0.17 


3.94±0.13 


3.79±0.16 


xVdof 


3266.3/3540 


4270.6/4298 


4284.8/4298 



also be reduced in strength. To understand the magni- 
tude of this effect, the Si+Au^ 2tt^ data sets were fit 
with a half-strength Coulomb correction, and also with 
no Co ulom b correction. The results are summarized in 
Table VII . Variations in the radii are as large as 20% for 
R\ and Ro, and A is greatly affected. A thorough investi- 
gation of effect of long lived resonances on the Coulomb 
correction is left for further study. 



APPENDIX C: 



TABLES OF CORRELATION 
RADII 



Tables VIII through XV show the complete set of fit- 
ted radii used to study the system, centrality, and trans- 
verse mass dependence for the Qinv, Yano-Koonin, and 
Bertsch-Pratt fits of the 14.6 A-GeV/c Si+Al, Si-^Au, 
and 11.6 A-GeV/c Au-I-Au data sets. All R and r pa- 
rameters are given in units of fm and fm/c, respectively. 
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TABLE VIII: Qinv fit parameters and < rriT > for data sets 1-9. 



System 


Set 


< niT > 


< > 


A 




xVdof 




1 
2 


0.348 
0.349 


0.327 
0.319 


0.491 ± 0.029 
0.676 ± 0.040 


3.59 ± 0.16 
4.42 ± 0.16 


32.1/27 
20.8/27 




3 


0.345 


0.323 


0.487 ± 0.039 


3.70 ± 0.23 


46.6/27 


Si+Au 


4 


0.345 


0.313 


0.551 ± 0.035 


4.47 ± 0.18 


45.0/27 


5 


0.342 


0.306 


0.511 ± 0.026 


4.91 ± 0.15 


39.2/27 




6 


0.349 


0.316 


0.514 ± 0.041 


4.54 ± 0.22 


11.5/27 




7 


0.411 


0.361 


0.441 ± 0.039 


5.13 ± 0.28 


24.1/27 


Au+Au 


8 


0.410 


0.348 


0.451 ± 0.051 


6.55 ± 0.43 


31.3/27 




9 


0.420 


0.358 


0.428 ± 0.046 


6.09 ± 0.38 


24.8/26 



TABLE IX; Yano-Koonin fit parameters for data sets 1-9. 



System 


Set 


Ym 


A 


R 


r 


P 


xVdof 


Si+Al 


1 


1.72 


0.452 + 0.024 


1.74 ± 0.13 


2.17 + 0.19 


-0.06 ± 0.06 


3770.7/4044 


2 


1.72 


0.627 + 0.035 


2.75 ± 0.14 


2.01 + 0.38 


-0.08 ± 0.07 


3901.9/3992 




3 


1.72 


0.467 + 0.024 


2.17 ± 0.11 


2.66 ± 0.27 


-0.04 ± 0.05 


3357.6/3605 


Si+Au 


4 
5 


1.72 
1.50 


0.503 + 0.028 
0.499 + 0.022 


2.45 ± 0.16 
2.83 ± 0.09 


2.70 ± 0.31 
3.01 ± 0.25 


-0.03 ± 0.06 
-0.06 ± 0.05 


4048.3/4156 
4588.6/4844 




6 


1.50 


0.515 + 0.038 


2.76 ± 0.20 


3.10 ± 0.39 


0.07 + 0.07 


3811.3/4055 




7 


1.45 


0.496 + 0.031 


2.84 ± 0.14 


3.86 ± 0.30 


0.01 ± 0.07 


4234.2/4407 


Au+Au 


8 


1.45 


0.454 ± 0.044 


3.36 ± 0.27 


4.48 ± 0.54 


-0.12 ± 0.08 


4382.3/4503 




9 


1.45 


0.530 ± 0.053 


3.96 ± 0.25 


4.23 ± 0.88 


0.11 ± 0.10 


4136.8/4486 



TABLE X; Bertsch-Pratt fit parameters for data sets 1- 
used in error propagation. 



■9. Values for r are calculated from Eq. 4.5 with Rs-Ro covariances 



System 


Sot 


A 


Ri 




Ro 


T 


xVdof 


Si+Al 


1 


0.437 ± 0.024 


1.76 ± 0.17 


1.05 ± 0.31 


2.61 ± 0.13 


2.71+0.14 


3507.3/3594 


2 


0.645 ± 0.045 


2.96 ± 0.23 


2.30 ± 0.24 


3.34 ± 0.17 


3.13+0.26 


3266.3/3540 




3 


0.452 ± 0.034 


2.25 ± 0.22 


1.63 ± 0.43 


3.21 ± 0.18 


3.06+0.10 


3053.2/3184 


Si+Au 


4 

5 


0.488 ± 0.036 
0.503 ± 0.021 


2.29 ± 0.21 
2.82 ± 0.14 


2.12 ± 0.30 
2.65 ± 0.18 


3.43 ± 0.18 
3.89 ± 0.13 


2.75+0.20 
3.32+0.20 


3382.1/3636 
4274.9/4298 




6 


0.517 ± 0.054 


2.57 ± 0.27 


2.56 ± 0.39 


3.97 ± 0.28 


3.52+0.31 


3544.5/3672 




7 


0.506 ± 0.037 


2.62 ± 0.23 


3.00 ± 0.33 


4.66 ± 0.26 


3.95+0.37 


2942.1/3114 


Au+Au 


8 


0.480 ± 0.043 


3.25 ± 0.32 


3.84 ± 0.46 


5.35 ± 0.33 


4.13+0.56 


3070.2/3236 




9 


0.536 ± 0.047 


3.92 ± 0.37 


3.79 ± 0.45 


5.60 ± 0.43 


4.55+0.61 


2694.4/3095 



TABLE XI: Qinv and Bertsch-Pratt fit parameters with _Rio cross-term (see Eq. 4 



for data sets 1-9. 



System 


Set 


A 


Ri 


Rs 


Ro 


Rio 


xVdof 


Si+Al 


1 


0.430 + 0.023 


1.65 + 0.17 


0.97 + 0.31 


2.59 + 0.13 


0.58 + 0.36 


3441.8/3604 


2 


0.618 + 0.034 


2.71 + 0.19 


2.20 + 0.23 


3.29 + 0.15 


0.51 + 0.69 


3363.7/3533 




3 


0.426 + 0.034 


2.03 + 0.23 


1.26 + 0.49 


3.15 + 0.24 


0.74 + 0.82 


3260.6/3200 


Si+Au 


4 


0.479 + 0.028 


2.16 + 0.18 


2.01 + 0.26 


3.41 + 0.16 


0.68 + 0.68 


3303.2/3643 


5 


0.496 + 0.021 


2.75 + 0.14 


2.59 + 0.18 


3.86 + 0.13 


2.59 + 0.75 


4126.7/4293 




6 


0.496 + 0.037 


2.47 + 0.24 


2.47 + 0.31 


3.93 + 0.24 


0.50 + 0.94 


3498.9/3621 




7 


0.499 + 0.036 


2.55 + 0.23 


2.90 + 0.32 


4.66 + 0.26 


-1.34 + 1.29 


2919.7/3170 


Au+Au 


8 


0.464 + 0.040 


3.18 + 0.37 


3.50 + 0.45 


5.37 + 0.22 


4.58 + 2.37 


3018.2/3264 




9 


0.582 + 5.033 


4.49 + 1.15 


3.71 + 1.00 


5.74 + 1.09 


-3.10 + 3.08 


3018.2/3264 
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TABLE XII: LCMS frame Bertsch-Pratt fit parameters for the for data sets 1-9. 



System 


Set 


A 






Ro 


2/1 r 

X /dof 


Si+Al 


1 




1 49 + n 99 








2 


0.53i ± 0.0 li 


2.()8 ± 0.25 


2.00 ± 0.27 


3.:i.5 ± 0.19 


:ii():i.2/:!i78 




3 


0.359 ± 0.038 


1.96 ± 0.30 


0.90 ± 0.82 


3.12 ± 0.26 


3113.5/3107 


Si+Au 


4 


0.392 ± 0.028 


1.97 ± 0.22 


1.57 ± 0.37 


3.41 ± 0.17 


3206.7/3551 


5 


0.398 ± 0.023 


2.66 ± 0.19 


2.23 ± 0.23 


4.05 ± 0.17 


3835.7/4199 




6 


0.413 ± 0.042 


2.43 ± 0.31 


2.11 ± 0.46 


4.17 ± 0.27 


3261.9/3447 




7 


0.399 ± 0.043 


2.28 ± 0.30 


2.54 ± 0.46 


4.77 ± 0.38 


3002.2/3084 


Au+Au 


8 


0.347 ± 0.037 


2.75 ± 0.36 


3.11 ± 0.57 


5.49 ± 0.39 


3034.3/3188 




9 


0.392 ± 0.042 


3.49 ± 0.44 


3.05 ± 0.56 


6.05 ± 0.58 


2714.5/3033 



TABLE XIII: Qinv and Bertsch-Pratt fit parameters for rriT-dependent data sets a-h. 



System 


Set 




< my > 


A 


Ri 


Rs 


Ro 


r 


xVdof 




a 


5.03 ± 0.21 


0.248 


0.560 ± 0.035 


3.52 ± 0.25 


3.25 ± 0.28 


4.33 ± 0.22 


3.32+0.35 


1834.6/1844 




b 


4.69 ± 0.29 


0.266 


0.577 ± 0.059 


3.10 ± 0.31 


3.15 ± 0.45 


4.01 ± 0.30 


2.75+0.54 


2036.9/2083 


Si+Au 


c 


4.84 ± 0.29 


0.333 


0.530 ± 0.064 


2.81 ± 0.29 


2.69 ± 0.34 


3.53 ± 0.25 


2.49+0.40 


2796.6/3007 




d 


3.94 ± 0.34 


0.408 


0.512 ± 0.075 


2.11 ± 0.28 


2.18 ± 0.44 


3.84 ± 0.46 


3.32+0.53 


3705.3/3657 




c 


3.93 ± 0.25 


0.441 


0.628 ± 0.058 


2.30 ± 0.17 


2.54 ± 0.25 


3.05 ± 0.20 


1.78+0.49 


4296.6/4230 




f 


6.05 ± 0.38 


0.294 


0.544 ± 0.055 


3.62 ± 0.38 


4.90 ± 0.51 


5.57 ± 0.41 


2.92+0.93 


1261.7/1215 


Au+Au 


g 


6.68 ± 0.85 


0.400 


0.479 ± 0.073 


3.05 ± 0.42 


3.11 ± 0.54 


4.39 ± 0.44 


3.32+0.64 


1960.6/2010 




h 


5.76 ± 1.09 


0.539 


0.508 ± 0.101 


2.36 ± 0.34 


3.73 ± 0.57 


4.21 ± 0.54 


2.10+1.44 


3541.3/3693 



TABLE XI\': Qi^ and Boitst-h-Pratt fit paraiuntnrs for low-in i data soth i 9. 



System 


Set 


< my > 


^inv 


A 


Ri 




Ro 


xVdof 




Si+Al 


1 

2 


0.277 
0.271 


3.38 + 0.19 
4.37 ± 0.20 


0.502 + 0.038 
0.685 + 0.044 


2.34 + 0.25 
3.36 + 0.25 


2.15 + 0.32 
2.58 + 0.28 


2.75 + 0.18 
3.46 + 0.20 


1867.9/1849 
1789.6/1854 




Si+Au 


3 
4 
5 
6 


0.270 
0.269 
0.264 
0.268 


4.23 + 0.31 
4.45 + 0.21 
5.01 + 0.17 
4.72 + 0.29 


0.500 + 0.070 
0.569 + 0.042 
0.539 + 0.027 
0.576 + 0.046 


2.59 + 0.45 
2.90 + 0.27 
3.02 + 0.18 
3.07 + 0.28 


2.30 + 0.63 
2.86 + 0.33 
3.18 + 0.22 
3.08 + 0.38 


3.39 + 0.36 
3.49 + 0.22 
4.05 + 0.17 
4.03 + 0.26 


1589.0/1699 
1796.6/1859 
2939.3/2981 
2105.1/2151 




Au+Au 


7 
8+9 


0.309 
0.308 


5.20 ± 0.33 
6.08 ± 0.37 


0.522 + 0.051 
0.542 + 0.041 


2.87 + 0.32 
3.61 + 0.30 


3.36 + 0.50 
4.55 + 0.42 


4.52 + 0.35 
5.41 + 0.33 


1305.1/1427 
1640.1/1532 




TABLE XV: Qinv and Bertsch-Pratt fit parameters for high-mT data sets 1-9. 


System 


Set 


< niT > 


Rinv 


A 


Ri 




Ro 


xVdof 




Si+Al 


1 
2 


0.410 
0.411 


3.58 ± 0.26 
4.28 ± 0.24 


0.429 + 0.039 
0.760 + 0.069 


1.37 + 0.41 
2.65 + 0.25 


0.13 + 0.78 
2.41 + 0.26 


2.31 + 0.35 
3.08 + 0.21 


3491.7/3590 
3428.1/3534 




Si+Au 


3 
4 
5 

6 


0.411 
0.411 
0.411 

0.412 


2.98 ± 0.23 
3.93 + 0.31 
4.31 ± 0.24 

3.90 + 0.34 


0.472 + 0.072 
0.519 + 0.052 
0.573 + 0.064 

0.495 + 0.079 


2.26 + 0.27 
1.97 + 0.19 
2.15 + 0.20 

2.14 + 0.26 


1.38 + 0.57 
1.82 + 0.31 
2.36 + 0.27 

2.14 + 0.49 


2.97 + 0.27 
3.17 + 0.19 
3.15 + 0.21 

3.72 + 0.42 


3187.4/3165 
3523.8/3633 
4840.6/4986 

3669.6/3654 




Au+Au 


7 
8+9 


0.485 
0.485 


4.24 + 0.47 
6.44 ± 0.85 


0.644 + 0.089 
0.465 + 0.056 


2.68 + 0.34 
2.34 + 0.28 


3.02 + 0.40 
3.41 + 0.39 


4.72 + 0.52 
4.01 + 0.35 


3174.5/3113 
3581.9/3701 
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